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airplane,  and  symmetry  can  be  assumed  for  the  other  half.  If, 
however,  the  free  stream  is  not  parallel  to  the  plane  of  symmetry 
and/or  if  an  instrument  pod  is  added  under  one  wing,  then 
symmetry  is  spoiled.  In  vortex  lattice  computations,  the  symmetry 
of  the  airplane  can  still  be  used  to  advantage,  not  only  in 
computation  of  the  coefficients  of  the  linear  algebraic  system, 
but  more  so  in  its  solution.  Computation  time  for  a  completely 
symmetric  system  is  reduced,  in  comparison  to  that  for  an 
asymmetric  system  of  the  same  size,  by  a  factor  of  eight  for 
one  plane  of  symmetry.  If  there  are  two  or  three  planes  of 
symmetry,  the  factors  are  64  and  512,  respectively.  If  the 
system  is  only  partially  symmetric,  the  corresponding  factors 
may  still  be  as  high  as  4 ,  16,  and  64. 
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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  of  the  research 
presented  were  obtained  by  ARO,  Inc.,  AEDC  Group  (a  Sverdrup  Corporation  Company), 
operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under 
ARO  Project  Number  P43T-07.  The  manuscript  w'as  submitted  for  publication  on 
November  7,  1980. 
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1.0  INTRODUCTION 

The  potential  flow  computer  program  (PFP)  described  in  Ref.  I  has  been  used 
frequently  for  aerodynamic  analyses  in  support  of  testing  at  the  Arnold  Engineering 
Development  Center  (AEDC).  In  one  particular  application,  the  investigation  of  the  gun 
gas  ingestion  for  the  A- 10,  the  airplane  was  symmetric  but  the  gun  was  off  center.  The 
completely  symmetric  model  required  four  hours  of  IBM  370/165  computer  time  to  run 
on  the  PFP.  To  obtain  a  solution  with  the  gun  off  center  would  have  required  about 
thirty  hours  per  run.  This  rush  project  motivated  a  quick  analysis  of  symmetry  properties 
which  resulted  in  a  procedure  whereby  the  off-center  runs  could  be  made  in  seven  hours. 
About  ten  such  runs  were  made.  After  the  project  was  completed,  a  general  method  was 
developed  as  a  spin-off  under  a  computational  fluid  dynamics  project. 

Symmetry  is  commonly  used  to  cut  computation  time;  however,  it  is  usually  an 
all-or-nothing  proposition.  This  need  not  be  the  case  for  vortex  lattice  computations.  If 
the  model  has  a  nonsymmetnc  portion  and/or  if  the  boundary  conditions  are  not 
symmetric,  the  symmetry  of  the  system  can  still  be  used  to  advantage,  not  only  in 
computation  of  the  coefficients  of  the  linear  algebraic  system,  but  also  in  its  solution. 
Computation  Time  for  a  completely  symmetric  system  is  reduced,  compared  to  that  for 
an  asymmetric  system  of  the  same  size,  by  a  factor  of  eight  for  one  plane  of  symmetry. 
If  there  are  two  or  three  planes  of  symmetry,  the  factors  are  64  and  512,  respectively.  It 
the  system  is  only  partially  symmetric,  the  corresponding  factors  may  still  be  as  high  as 
4.  16.  and  64. 

Section  2.0  presents  the  method  for  taking  advantage  of  partial  symmetry.  Beginning 
with  the  pertinent  information  on  the  vortex  lattice  method,  it  shows  how  the  algebraic 
system  representing  a  symmetric  geometry  satisfies  certain  symmetry  conditions.  If  the 
geometry  is  only  partially  symmetric,  the  algebraic  system  can  be  partitioned  so  that  one 
portion  satishes  the  symmetry  conditions.  The  symmetry  conditions  can  be  used  to 
produce  a  quicker  solution  of  the  symmetric  portion  of  the  system.  Section  3.0 
summarizes  and  does  an  operation  count  which  verifies  the  value  ot  the  method. 

In  the  present  domain  of  application,  the  number  ot  planes  of  symmetry,  X,  is 
meaningful  only  for  X  =  1.  2.  or  3.  Mathematically,  however,  the  method  generalizes  to  X 
being  an  arbitrary  positive  integer.  The  validity  of  the  method  and  existence  theorems  on 
jhe  inverses  involved  can  then  be  proved  by  mathematical  induction  This  mathematical 
analysis  is  presented  in  Lite  Appendix. 
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2.0  METHOD 


2.1  VORTEX  LATTICE  METHOD 

The  vortex  lattice  method  (Ref.  1)  models  the  boundary  of  a  flow  regime  with  N 
singularities,  either  vortices  or  sources.  Each  singularity  induces  flow  with  velocities 
proportional  to  its  strength.  The  resultant  flow  is  the  superposition  of  the  flows 
corresponding  to  all  the  singularities  and  an  optional  free  stream.  The  strengths  of  the  N 
singularities  are  determined  by  imposing  boundary  conditions  at  N  points  called  control 
points.  This  produces  a  linear  system  of  N  algebraic  equations 

A  X  =  R  (1) 

to  be  solved  for  the  unknown  strengths,  X. 

The  equations  for  the  coefficient  matrix,  A,  and  the  right-hand  side,  R,  are  not 
pertinent  to  this  analysis,  but  the  fact  that  the  coefficient  matrix  is  determined  only  by 
the  geometry  of  the  boundary  is  crucial  to  the  development.  This  fact  will  be  used  in 
Section  2.3  to  derive  symmetry  conditions  which  are  present  in  A  if  symmetry  exists. 

Assume  that  the  boundary  has  a  geometrically  symmetric  portion  made  up  of  M 
singularities  and  an  asymmetric  portion  made  up  of  p  singularities,  so 

M  =  M  +  p  (2) 

then  the  system  of  equations,  Eq.  (1),  can  be  partitioned  to  take  the  form 


(3) 


where  A  is  M  x  M,  B  is  M  x  p,  C  is  p  x  M,  D  is  p  x  p,  X  is  M  x  1 ,  Y  is  p  x  1,  R  is  M  x 
1,  and  S  is  p  x  1.  Under  such  a  partition,  the  matrix  A  can  be  constructed  so  that  it 
satisfies  the  symmetry  conditions. 

Equation  (3)  can  be  written 


AX  +  BY  =  R  (4) 

and 

CX  +  DY  =  S  (5) 
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The  purpose  of  the  method  being  presented  is  to  use  the  symmetry  conditions  present  in 
A  to  effect  a  fast  solution  of  Eq.  (4)  for  X  in  terms  of  Y.  obtaining 

X  =  r  -  bY  (6) 

Equation  (6)  can  then  be  substituted  into  Eq.  (5)  to  obtain 

(D  -  Cb)  Y  =  IS  —  Cr)  (7) 

and  Gaussian  elimination  can  be  used  on  Eq.  (7)  to  solve  for  Y.  The  Y  can  be  substituted 
into  Eq,  (6)  for  the  complete  solution. 

Some  discussion  is  needed  on  the  existence  of  the  solutions.  Equation  (1)  has  a 
solution  if  A  is  nonsingular.  The  step  from  Eq.  (4)  to  Eq.  (6)  can  be  taken  if  A  is 
nonsingular.  If  A  and  A  are  nonsingular,  then  the  matrix,  (D  —  Cb),  of  Eq.  (7)  is 
nonsingular.  The  matrix,  A.  will  be  nonsingular  if  the  system  has  been  modeled  correctly. 
The  matrix.  A.  will  be  nonsingular  if  the  symmetric  portion  has  been  modeled  correctly. 
In  general,  if  a  coefficient  matrix  is  singular  or  even  ill-conditioned  it  is  because  the 
system  has  been  improperly  modeled. 

2.2  NUMBERING  THE  SINGULARITIES 

The  vortex  mode!  is  composed  of  either  vortices  or  sources  commonly  called 
singularities  because  the  corresponding  velocity  functions  are  singular  at  their  centers  (not 
to  be  confused  with  matrix  singularity).  The  partitioning  of  the  system  depends  on  the 
ordering  of  these  singularities.  The  following  numbering  scheme  will  be  used-  Let  A  be 
the  number  of  planes  of  symmetry,  A  =  1 ,  2,  or  3.  It  A  =  1 ,  then  the  symmetric  portion 
of  the  vortex  model  can  be  divided  into  two  parts.  One  part  will  be  called  the  basis  and 
the  other  its  reflection.  If  A  =  2,  then  the  symmetric  portion  can  be  divided  into  four 
parts,  one  basis  and  three  reflections.  If  A  =  3,  there  are  eight  parts  with  one  basis  and 
seven  reflections.  Let  A  be  the  number  of  parts,  then 

A  =  2*-  (8) 

Let  n  be  the  number  of  singularities  in  the  basis:  then  the  number  of  singularities  in  the 
symmetric  portion  is 


M  =  An  (.9) 

Let  the  planes  of  symmetry  be  represented  by  Pj.  where  H  =  1  co  A.  The  singularities  ol 
the  basis  will  be  numbered  from  l  to  n.  The  reflection  ot  the  i  111  singularity.  1  <  i  <  n. 
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with  respect  to  PL  will  be  numbered  n  +  l.  For  X  >  L  the  reflection  of  the  i'tli 
singularity,  1  <  i  <  2n,  with  respect  to  P2  will  be  numbered  2n  +  i.  For  X  >  2,  the 
reflection  of  the  i  th  singularity,  1  ^  1  ^  4n  with  respect  to  P3  will  be  numbered  4n  +  1. 
As  already  implied  by  the  form  of  Eq.  (3),  the  nonsymmetric  portion  will  be  numbered 
from  M  +  I  to  N. 

With  this  numbering  system,  the  matiices,  A.  X,  B.  and  R  can  be  partitioned  with 
respect  to  the  A  parts 

A  =  ((A,j))  X  =  ((X,» 

B  =  {(B,.))  R  =  ((R,))  U0) 


where  i.  J  -  1  to  A.  A^  is  n  x  n,  X,  is  11  x  3 ,  B,  is  11  x  p,  and  R,  is  n  x  1.  With  this 
partition,  Eq.  (4)  can  be  written 


A 

I,  A,j  X;  +  B,  Y  =  R,  (11) 


with  i  -  I  to  A.  As  a  special  case,  if  there  is  no  nonsymmetric  portion,  then  Eq.  (11) 
takes  the  form 

A 

]f]  Au  xi  =  R>  (12) 

with  i  -  1  to  A  and  there  will  be  no  Eq.  (5).  This  report  offers  no  new  approach  to  the 
other  special  case  of  no  symmetric  portion,  therefore,  it  will  not  be  considered. 

Regarding  notation,  it  will  at  times  be  convenient  to  add  a  superscript  to  a  variable 
indicating  the  number  of  planes  of  symmetry,  X.  This  appendage  will  m  no  way  alter  the 
definition  of  the  variable. 

2.3  SYMMETRY  CONDITIONS  FOR  A 

Since  A  is  the  coefficient  matrix  for  the  symmetric  portion  and  it  depends  only  on 
the  boundary  of  the  system,  certain  deductions  about  its  internal  form  can  be  made. 
These  symmetry  conditions  will  now  be  derived. 

Each  part  of  the  symmetric  portion  has  a  corresponding  part  with  respect  to  each 
plane,  the  corresponding  part  being  its  reflection  with  respect  to  that  plane.  A  table  can 
be  made  of  these  corresponding  parts,  which  is  determined  by  the  numbering  system  of 
Sec.  2.2. 
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Table  1.  Corresponding  Reflected  Parts 


Plane  A  =  1  A  =  2  A  =  3 


1 

3 

4 

5 

6 

7 

8 

Pi 

2 

1 

4 

3 

6 

5 

8 

7 

P  2 

3 

4 

1 

7 

8 

5 

6 

Pi 

5 

6 

7 

8 

1 

O 

jL 

3 

4 

The  table  for  A  =  1  is  above  and  to  the  left  of  the  first  bold  line,  and  the  table  for 
A  =  2  is  above  and  to  the  left  of  the  second  bold  line.  The  third  bold  line  marks  the 
table  for  A  =  3. 


A  given  matrix.  A(J .  depends  only  on  the  geometry  of  parts  i  and  j  and  their  relative 
position.  Given  a  plane  of  symmetry,  Pg,  if  i'  and  j'  are  the  corresponding  parts  of  i  and  j 
with  respect  to  Pg,  then  parts  i'  and  j'  (since  they  are  reflections)  have  the  same 
geometry  and  relative  positions  as  do  parts  i  and  j.  Therefore,  A,1/  =  AtJ.  Such  a  relation 
can  be  written  lor  each  submatrix,  AIJ(  of  A  with  respect  to  each  plane  of  symmetry,  Pg. 
Thus,  there  aTe  A(Aj2  such  relations.  For  example,  when  A  =  2,  the  relations  for  the  first 
row  of  A  are 

Pi  Relations  P;  Relations 

All  =  At  2  An  =  A3  3 


A] 2  =  At  1 
A] 3  =  A24 
A] 4  =  A; 3 


A  l  2  =  A3  4 
Al  3  “  A3  1 

A14  =  A3  2 


For  A  =  2,  there  are  24  more  relations  which  have  not  been  written  out.  For  a  given  A,  if 
all  the  relations  (not  all  are  independent)  are  written  out,  A  can  be  written  in  terms  of 
the  submatrices,  A^,  of  the  first  row.  Define 


Aj  -  Aij 

Then  the  form  of  A  for  A  -  1,2,  and  3  is 


(13) 
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Ai  A2 

A1  =  (14) 

A2  A] 


A  pattern  can  be  noticed  in  Ak;  namely,  if  A2  and  A3  are  partitioned  as  shown,  then 
each  partition,  large  or  small,  has  the  same  form  as  the  right-hand  side  of  Eq.  (14).  This 
is  the  essence  of  a  recursion  procedure  for  generating  A*  in  the  Appendix. 

2.4  SYMMETRY  CONDITIONS  FOR  R 

It  is  not  necessary  for  the  boundary  conditions  to  be  symmetric  to  use  the  solution 
of  Section  2.5.  In  fact,  symmetry  of  the  boundary  conditions  is  of  no  practical  advantage 
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except  in  tlic  special  case,  p  =  o.  In  that  case,  if  R  satisfies  the  symmetry  conditions, 
there  is  an  additional  timesaving  factor  of  A  in  the  solution. 

Boundary  conditions,  with  respect  to  a  plane,  can  be  either  symmetric  (S-sy  mine  try) 
or  antisymmetric  (A-sy  mine  try).  As  an  illustration,  for  X  —  1  il  a  plane  is  symmetric,  then 


R2  =  R| 

and  if  it  is  antisymmetric,  then 

R2  =  -Ri 


Since  each  plane  can  be  either  symmetric  or  antisymmetric,  there  are  A  possible 
combinations.  The  possibilities  are  shown  in  Table  2. 

Table  2,  Types  of  Symmetry 


'  Plane  X  =  1  X  =  2  X  -  3 


1 

Pe 

1 

-) 

3 

4 

5 

6 

7 

8 

Pi 

s 

A 

S 

A 

S 

A 

S 

A 

Pi 

s 

S 

A 

A 

S 

s 

■ 

A 

A 

P3 

s 

S 

S 

S 

A 

A 

A 

A 

The  convention  for  using  Table  2  for  X  =  1.2,  or  3  is  above  and  to  the  left  of  the 
indicated  bold  line,  and  is  the  same  as  explained  for  Table  1.  Table  2  can  be  used  to 
construct  a  reflection  matrix 


n  =  (o?u>)  (17) 

i.  j  =  1  to  A.  The  elements  of  77  are  plus  or  minus  ones.  The  first  row  consists  of  plus 
ones.  The  second  row  of  17  is  determined  from  the  first  row  and  the  P|  row  ot  Table  2. 
An  S  in  Table  2  means  copy  fiom  the  first  row.  and  an  A  means  copy  with  a  change  of 
sign.  This  process  is  analogous  to  the  reflection  of  the  basis  with  respect  to  l\ .  For  X  > 
1,  the  third  and  fourth  rows  of  17  arc  determined  from  the  first  two  rows  and  the  P?  row 
of  Table  2.  An  S  in  Table  2  means  copy  from  the  first  two  rows  and  an  A  means  copy 
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with  a  change  of  sign.  This  is  analogous  to  the  reflection  of  the  first  two  parts  with 
respect  to  P2  to  obtain  parts  three  and  four.  For  X  >  2,  rows  five  through  eight  are 
obtained  from  the  first  four  rows  and  the  P3  row  of  Table  2  in  like  manner:  this  is 
analogous  to  the  reflection  of  the  first  four  parts  to  obtain  parts  five  through  eight.  By 
these  rules,  the  following  is  obtained: 


V  = 


-l 


and 


1111 
1  -11-1 

1  1-1-1 


1 


-1  -1 


1 


1  1111 

1-11-11 
1  I  -1  -I  1 
1  -1-111 

1  111-1 

1  -11-1  -1 

1  1-1-1  -1 


! 


-1  -1 


I  -1 


1  1 

I  -1 

-1  -1 

-1  1 

-I  -1 


-1  1 

1  ! 

1  -1 


Some  useful  properties  of  77  are  noted.  It  is  symmetric 

17*  =  r) 


(18) 


(19) 


(20) 


(21) 
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where  the  asterisk  indicates  the  transpose.  Its  inverse  is  given  by 

i'1  =  x  " 


(22) 


This  is  proved  in  Theorem  1  of  the  Appendix.  The  reflection  matrix  is  an  example  of  a 
Hilda mard  matrix.  (Refs.  2,  3.  and  4).  A  Hadamard  matrix  is  a  square  matrix  of  plus  and 
minus  ones  such  that 

HH*  =  m  I 


where  m  is  the  order  of  the  matrix. 

Having  available  the  reflection  matrix,  the  symmetry  conditions  for  R  are  simply 
stated.  If  the  boundary  conditions  have  symmetry  of  the  j'th  type,  then 


R,  =  i7,j  Rt 


(23) 


for  i  =  I  to  A. 


2.5  SOLUTION 

As  stated  in  Section  2.1.  the  objective  of  this  report  is  a  quick  solution  of  Eq.  (4). 
Sections  2.2.  2.3.  and  2.4  were  in' preparation  for  this  solution.  To  this  end.  Eq.  (11)  is 
multiplied  by  T?kj  and  summed  with  respect  to  i  to  obtain 


a  \ 

£  £ 

.=  1  j=  L 


A 

T?ki  A, |  X,  +  £ 

1=  I 


Pki  ^  Z  nki  Ri 

i- 1 


for  k  =  1  to  A.  Define 


akj  X|  1  Aij 
=  “  j£  nkj  x, 

i  A 

/3k  ^  ~  Pkj  B, 


(24) 


(25) 


(26) 


(27) 


]  A 

Pk  “  ~  Z  Pk  i  B, 


(28) 
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and  the  corresponding  matrices 

a.=  ((akj))  6  =  ((0k)) 

0  =  ((Ac))  P  =  ((Pk)) 

with  k,  j  =  1  to  A. 

Fundamental  to  the  solution  is 


(29) 


akj  ~  Pkj  a-k  i  (30) 

k,  j  =  1  to  A.  This  can  be  verified  by  obtaining  the  product  t?A.  Except  for  the  sign,  the 
elements  of  a  given  row  are  the  same,  and  the  signs  of  corresponding  elements  of  ij  and  a 

are  the  same.  The  formal  proof  of  Eq.  (30)  is  given  in  Theorem  3  of  the  Appendix.  Using 

Eqs.  (25)  through  (30),  Eq.  (24)  becomes 

Okt  #k  +  Pk  Y  =  pk  (31) 

for  k  -  1  to  A.  Thus,  instead  of  solving  Eq.  (4)  which  has  M  rows,  there  are  A  systems 

of  the  form  of  Eq.  (3 1 ),  having  n  rows  each.  This  can  result  in  a  considerable  savings,  as 
will  be  seen  in  Section  3.0.  Equation  (31)  can  be  solved  for  0k  if  akj  is  nonsingular.  In 
Theorem  4  of  the  Appendix,  aki,  k  =  1  to  A  are  shown  to  be  nonsingular  if  A  is 
nonsingular.  Thus,  Qk  can  be  written  in  terms  of  Y 


9k  =  rk  -  bk  Y 


(32) 


for  k  -  1  to  A.  Numerically,  it  is  much  more  efficient  to  form  the  matrix 


<*k  1 


Pk 


j  Ac] 


(33) 


and  solve  by  Gaussian  elimination  than  to  invert  ak  i . 


Having  6  in  terms  of  Y,  X  can  be  obtained  in  terms  of  Y  from  Eqs.  (26)  and  (22): 


Xj  2  77jk  0k 
k=  i 


(34) 


i  1  to  A.  From  Eqs.  (10),  (34),  and  (32),  the  r  and  b  of  Eq.  (6)  can  be  assembled  as 

r  =^(J1  mk  7k|  (35) 
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and 

b  -((&  0..  W)) 

This  completes  the  solution  promised  in  Section  2.1. 


(36) 


There  remains  the  special  case,  p  =  o.  Applying  the  same  procedure  to  Eq.  (12) 
produces 


«k  i  Ok  =  pk  (37) 

instead  of  Eq  (31)  for  k  =  1  to  A  There  is  an  added  benefit  if  R  satisfies  the  symmetry 
conditions.  Assume  the  boundary  conditions  have  symmetry  of  thej'th  type,  substituting  Eq. 
(23)  into  (28),  and  using  Eq.  (22)  gives 


Pk  =  6kj  Ri  (38) 

k  =■  1  to  A.  where  5k j  is  the  Kronccker  Delta.  So  p|<  =  o,  implying  6^  =  o  for  all  k 
except  k  =  j.  Thus,  Eq.  (37)  needs  to  be  solved  only  for  k  =  j.  As  a  side  note,  the 
boundary  conditions  do  not  have  to  be  symmetric  with  respect  to  all  the  planes  of 
symmetry  to  be  beneficial.  Each  plane  for  which  the  boundary  conditions  have  either 
S-symmetry  or  A-symmetry  eliminates  half  of  the  systems,  Eq.  (37). 

3.0  SUMMARY 

The  algorithm  can  be  described  as  follows.  Given  A,  R,  X,  n,  and  p,  the  following 
steps  are  performed. 

1.  Compute  a.  Eq.  (25). 

2.  Compute  (3.  Eq.  (27). 

3.  Compute  p,  Eq.  (28). 

4.  For  k  =  1  to  A,  apply  Gaussian  elimination  to  the  matrix  of  Eq.  (33). 

5.  Compute  r,  Eq.  (35). 

6.  Compute  b,  Elq.  (36). 

7.  Compute  the  system  of  Eq.  (7). 
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8.  Solve  the  system  for  Y  using  Gaussian  elimination. 

9.  Compute  X,  Eq.  (6). 

The  time  required  to  perform  this  task  depends  on  the  number  of  operations,  i.e., 
additions,  subtractions,  multiplications,  and  divisions.  Different  operations  require 
different  times,  and  even  the  ratios  of  times  vary  with  computer.  Considering  parallel 
processors,  how  well  a  method  vectorizes  must  be  taken  into  account.  Such  a  detailed 
accounting  seemed  to  be  a  diversion,  so  the  rudimentary  convention  of  Ref.  5,  counting 
only  multiplications  and  divisions,  will  be  followed. 

The  first  savings  from  symmetry  come  in  the  computation  of  the  coefficient  matrix. 
Equations  (14)  through  (16)  show  that  only  A  of  the  (A)2  submatrices  of  A  need  to  be 
calculated.  Similarly,  if  R  satisfies  the  symmetry  conditions,  Eq.  (23),  only  one  of  the  A 
submatrices  of  R  need  be  calculated. 

In  the  solution.  Step  1  requires  only  additions  and  subtractions,  with  tj  just 
supplying  the  signs.  The  same  can  be  said  for  Steps  2  and  3,  except  Step  2  has  Mp 
divisions  by  A,  and  Step  3  has  M  divisions  by  A.  According  to  Ref.  5,  the  number  of 
operations  to  perform  Step  4  is 

A  n3  +  (p  +  1)  n2  —  y  n  J 

Steps  5  and  6  require  only  additions  and  subtractions.  Ih  Step  7,  the  computation  of  Cb 
requires  Mp2  multiplications,  and  the  computation  of  Cr  requires  Mp  multiplications. 
Step  8  requires 

y  P3  +  P2  -  y  P 

operations.  Step  9  requires  Mp  multiplications.  The  total  number  of  operations,  N, ,  is 

Ni  =  ~r  An3  +  An2  +  -f-  An 
+  Amp  +  Anp2  +  3Anp 

+  -j-  p3  +  p2  -  -y  p 

If  Eq.  (1)  were  solved  directly,  the  number  of  operations,  N2.  would  be 

l  1 

N?  =  —  (M  +  p)3  +  (M  +  p)2  +  —  (M  +  p) 

or 
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NN  =  — -  A3n3  +  A2n2  -  An 
3  3 

+  A2n2p  +  Anp2  +  2Anp  (40) 

+  -j-  p2  +  p-  -  —  p 

The  number  of  operations  saved  is 

N,  -  Nj  =  -j  (A  “  I)  A  (A  +  1)  n3 

+  An  [(A  -  1)  n  -  l]  (p  +  1)  (41) 

tqimtion  (41)  is  a  measure  of  the  value  of  exploiting  partial  symmetry. 

Often,  exploiting  symmetry  is  not  just  a  timesaver.  Many  times  it  is  a  matter  of 
whether  or  not  a  job  can  be  done;  the  factor  of  8*  pushes  it  beyond  the  realm  of 
practicality. 
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APPENDIX  A 

MATHEMATICAL  ANALYSIS 

Definition  1:  Super  Symmetric  Block  Matrix  (SSBM) 

Given  n  >  1 ,  X  o  and 

S  =  is  an  n  x  n  matrix} 

A  sequence  of  set 
procedure 

•  S0  =  S 


nifl  = 


,  j 

=  0 

to  X  can  be 

A 

A 

A  e  Sj.]j 

A 

A 

A  e  Sj.i  ( 

j  =  1  to  X.  If  AJeS^,  then  AJ  is  called  a  super  symmetric  block  matrix  of  order 
j.  It  can  be  noted  that 


•  If  Ak  and  Ak  are  SSBM's  of  order  k.  then 


Ak+1 


Ak 

Ak 

Ak 

Ak 

(A- 1 ) 


is  an  SSBM  of  order  k  +  1,  and  Ak  +  Ak  and  Ak  -  Ak  arc  SSBM'S  of  order  k. 
•  A1 .  A2,  and  A3  of  F.qs.  (14)  through  (16)  are  SSBM'S  of  orders  1.  2.  and  3, 
respectively. 


Definition  2.  The  Reflection  Matrix 

Given  X  >  o.  A  sequence  of  matrices  77J.  j  =  o  to  X,  can  be  defined  using  the 
following  recursive  procedure- 


r?°  =  1 


V' 


rjJ ' 1  17]  ‘ 1 

T|-1-1  -TfJ'1 


(A-2) 
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j  =  1  to  X.  The  matrix  77*  is  called  the  reflection  matrix  of  order  X.  It  can  be 
verified  that  rj1 ,  rj2 ,  and  p3  of  Eqs.  (18)  through  (20)  are  the  reflection 
matrices  of  order  1,  2,  and  3,  respectively.  It  is  easily  proved  that  the  first  row 
and  column  ol  a  reflection  matrix  consist  of  plus  ones. 

Lemma:  Let 


E  = 


F  F 
G  ~G 


(A-3) 


where  F  and  G  are  arbitrary  n  x  11  matrices.  The  block  matrix,  E,  is  nonsingular  if  and 
only  if  F  and  G  are  nonsingular;  if  nonsingular,  then 


Proof:  Note  that 


E  1 


F-i  C-l 

x 

2  F'1  -G-1 


F  F 
G  -G 


where  I  is  an  n  x  n  identity  matrix  and  0  is  an  n  x  n  zero  matrix. 


(A-4) 


Since 

I  I 

I  -1 

is  nonsingular,  E  is  nonsingular  if  and  only  if  F  and  G  are  nonsingular.  If  F  and  G  are 
nonsingular,  Eq,  (A-4)  can  be  easily  verified. 

Theorem  1 .  The  reflection  matrix  of  order  X  is  nonsingular  and  its  inverse  is  given  by 
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b 


where  A  =  2K. 

Proof'  The  theorem  is  obviously  true  for  X  =  o  since  A  —  r?  =  1  Assume  the  theorem  is 
true  for  X  =  k  for  some  k  >  o.  Applying  the  Lemma  to  Eq.  (A-2)  and  substituting  Eq. 
(A-5)  into  the  result  produces 


(r?fc  4  ’)' 1 


2-k  r?k  2'k  ^ 
2-k  rjk  -2-k  77k 


Factoring  out  the  2~k.  it  is  seen  that  the  theorem  is  true  for  X  =  k  ■+  1  and  thus  for  all  X. 

There  is  occasion  in  the  following  analysis  to  take  the  product  of  77  and  a  block 
matrix.  The  elements  of  one  are  scalars  and  the  elements  of  the  other  are  n  x  n  matrices. 
Since  the  product  of  a  scalar  and  a  matrix  is  defined,  the  product  of  tj  with  a  block 
matrix  can  be  defined  in  the  normal  manner. 

Theorem  2:  A  block  matrix.  A,  with  n  \  n  blocks  is  an  SSBM  of  order  X  if  and  only  if 


or  bv  Theorem  1 


A  =  (17^)  1  D  77^ 


A  =  —  77^  D 
A 


where  D  is  a  block  diagonal  matrix,  that  is 


D  = 


D, 


D: 

o''  D, 


(A-6) 


(A-7) 


(A-8) 


where  D, ,  j  =  1  to  A.  are  n  x  n. 

Proof  The  theorem  is  obviously  true  for  X  =  o,  since  ri°  =  1  so  A  =  D.  Assume  the 
theorem  is  true  for  X  =  k  for  some  k  >  o.  Given  an  arbitrary  D  of  the  form  of  Eq.  (A-8) 
with  A  =  2k  leads  to  the  definition 

D  =  Diag  (Dj ,  .  .  .  .  DttA)  ( A-9) 


and 
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D  =  Diag  (Di/,  A  +  1 .  .  .  .  ,  Da) 


(A- 10) 


and  Hq.  (A-8)  can  be  written 


0 

D 


Equation  (A-7)  then  becomes 


1 

T?k  7Jk 

D 

0 

TJk  7fk 

A  2k+1 

0 

5 

7?k  ~??k 

(A- 1  1) 


(A-l  2) 


which  can  be  wntten 


where 


and 


(A-13) 


(A-14) 


(A-l  5) 


The  quantities  in  parentheses  are  obviously  block  diagonal  matrices;  thus,  since  the 
theorem  is  assumed  true  for  X  =  k,  the  matrices  A  and  A  are  SSBM's  of  order  k.  Thus, 
by  the  form  of  Eq.  (A-13),  the  matrix  A  is  an  SSBM  of  order  X  ~  k  +  I. 


Assume  A  is  an  SSBM  of  order  X  -  k  +  I ;  then  A  can  be  written 


A 

A 


where  A  and  A  are  SSBM's  of  order  X  -  k.  Since  the  theorem  is  assumed  to  be  true  for  X 

=  K 
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and 


A  =  7?k  D  T?k 


A  =  —  J?k  D  Tjk 


where  D  and  D  are  block  diagonal  matrices.  By  matrix  algebra,  it  can  be 

i 


A  - 


?k+  1 


V 


k+1  D  1?k+l 


where 


D  +  D  0 
D  "  0  D-D 

Thus,  assuming  the  truth  of  the  theorem  for  X  =  k,  it  is  true  for  X  =•  k 
all  X. 

Corollary  if  A  is  a  nonsingular  SSBM,  then  so  is  its  inverse. 

Proof:  Since  A  is  an  SSBM 

A  =  TJ D  7} 

where  D  is  a  block  diagonal  matrix.  Its  inverse  is 

A-1  =  rj'1  D_l  t? 

Since  D’1  is  also  a  block  diagonal  matrix,  A-1  is  an  SSBM. 

Theorem  3’  Given 

a  =  (fcLjk ))  =  t?A 


with  j,  k  =  I  to  A  where  A  is  an  SSBM.  Then 

djk  =  hjk  djl 

Proof-  Since  A  is  an  SSBM.  substituting  Eq.  (A-6)  into  (A- 16)  produces 

a  =  77617' 1  D  ?})  =  Dij 


verified  that 


+  1 ,  and  thus  for 


(A- 16) 


(A-l  7) 


O 
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where  D  -  Diag  (Dj  ,....,DA )  with  Dj,  j  -  1  to  A,  being  n  x  n  matrices.  Therefore, 


ajk  -  (D*?)jk  =  Dj  77,  k 


(A- 1 8) 


when  k  =  1 


aji  =  Dj  =  dj  (A-19) 

since  the  first  column  of  i?  consists  of  plus  ones.  Subsituting  Eq.  (A-19)  into  (A-18) 
produces  the  desired  result,  Eq.  (A- 1 7). 

Theorem  4:  Given 


a  =  rjA  (A-20) 

where  A  is  a  nonsingular  SSBM,  the  block  matrix  a  is  obviously  nonsingular,  so  let 

■  x  «*»  (A-2» 

Then, 


a)k  =  ak’j 

or  by  Theorem  3 


(A-22) 


ajk  -  V)k  a-k’t 


(A-23) 


Proof1  Since  A  is  an  SSBM, 


a  =  17(77- 1  D  17)  =  D  i? 
and 

o'1  =t?'*  D1  =  —  17  D-1 
A  ' 

Therefore,  by  the  form  of  Eq.  (A-21) 


Qjk  =  (77  D1 1  )j k 


or 


°jk 


(A-24) 


The  desired  result.  Eq.  (A-23)  is  created  when  Eq.  (A-19)  is  substituted  mto  Eq.  (A-24). 


# 
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NOMENCLATURE 

A  Coefficient  matrix  for  the  symmetric  portion  of  the  model.  Eq.  (3) 

A,j  Submatrix  of  A,  Eq.  (10) 

A,  Atj,Eq.  (13) 

A  Coefficient  matrix  for  the  whole  system,  Eq.  (1) 

B  Submatnx  of  A,  Eq.  (3) 

B,  Submatrix  of  B,  Eq.  (10) 

b  Eqs.  (6)  and  (36) 

5k  Eq-  (32) 

C  Submatrix  of  A,  Eq.  (3) 

D  Coefficient  matrix  for  the  nonsymmetric  portion  of  the  model,  Eq.  (3).  In  the 

Appendix,  D  is  a  block  diagonal  matrix 

I  Identity  matrix 

M  Number  of  singularities  in  the  symmetric  portion  of  the  model,  Eq.  (9) 

N  Number  of  singularities  in  the  whole  model,  Eq.  (2) 

N]  Operation  count  using  symmetry,  Eq.  (39) 

Nj  Operation  count  neglecting  symmetry,  Eq.  (40) 

n  Number  of  singularities  in  the  basis  of  the  symmetric  portion  of  the  model 

Pg  C'tli  plane  of  symmetry 

p  Number  of  singularities  in  the  nonsymmetric  portion  of  the  model 

R  Right-hand  side  for  the  whole  system,  Eq.  (1) 

r  Eqs.  (6)  and  (35) 

R,  Right-hand  side  for  Part  i,  Eq.  (10) 

rk  Eq.  (32) 

R  Right-hand  side  for  the  symmetric  portion  of  the  model,  Eq.  (3) 

S  Right-hand  side  for  the  nonsymmetric  portion  of  the  model,  Eq.  (3) 

X  Strengths  of  the  symmetric  portion  of  the  model,  Eq.  (3) 
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X,  Strengths  of  Part  i,  Eq.  (10) 

X  Strengths  of  the  whole  system,  Eq.  (I) 

Y  Strengths  of  the  nonsymmotrie  portion  of  the  model,  Eq.  (3) 

a  Eq.  (29) 

akj  Eq.  (25)  ' 

P  Eq.  (29) 

Pk  Eq.  (27) 

8„  Kronecker  Delta,  Eq.  (38) 

r?  Reflection  matrix,  Section  2.4 

rj|j  Elements  of  tj,  Eq.  (17) 

6  Eq. (29) 

0k  Eq. (26) 

A  Number  of  parts,  Eq.  (8) 

X  Number  of  planes  of  symmetry 

P  Eq.  (29) 

Pk  Eq.  (28) 

NOTATION 

((  ))  Denotes  a  matrix  with  the  elements  indicated 

SYMBOLS 

*  Denotes  the  transpose 

-1  Denotes  the  inverse 

Any  superscript  on  a  variable  other  than  V'  or  ”-l"  does  not  alter  the  meaning  of  the 
variable,  but  refers  to  the  number  of  planes  of  symmetry  ,X. 


26 


